The Origin of the Decoupling of Oxygen and Silicon Dynamics in Liquid Silica as 
Expressed by its Potential Energy Landscape 
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The oxygen and silicon dynamics in silica is compared via computer simulations. In agreement 
with experimental data and previous simulations a decoupling of oxygen and silicon dynamics is 
observed upon cooling. The origin of this decoupling is studied in the framework of the potential 
energy landscape. From analysis of the transition features between neighboring superstructures 
of minima, denoted metabasins, the differences between the oxygen and the silicon dynamics can 
be quantified. The decoupling can be explicitly related to the presence of generalized rotational 
processes, giving rise to oxygen but not to silicon displacement. Closer analysis of these processes 
yields important insight into the nature of the potential energy landscape of silica. The physical 
picture of relaxation processes in silica, obtained in previous work for the oxygen dynamics, is 
consistent with the decoupling effects, elucidated here. 

PACS numbers: 64.70.Pf,65.40.Gr,66.20.+d 
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I. INTRODUCTION 

Silica is a well-known prototype network glass former. 
It is classified as a strong liquid in the Angell scheme p) 
because the temperature dependence of its dynamics dis- 
plays an Arrhcnius behavior. Interestingly, oxygen and 
silicon dynamics are characterized by different activation 
energies. Close to the glass transition temperature T g the 
difference between the self-diffusion constants for oxygen 
Do and silicon D$i is nearly two orders of magnitude 
0,0. Moreover, it has been shown that the viscous dy- 
namics of liquid silica close to T g is strongly controlled by 
the silicon dynamics Q. The large difference in diffusiv- 
ity indicates different relaxation mechanisms for oxygen 
and silicon atoms in liquid silica at low temperatures. 

The decoupling of silicon and oxygen single-particle 
dynamics has been already studied via computer simula- 
tions of Si0 2 @ modelled by the BKS potential [(J. It 
turned out that in the limit of very high temperatures 
the diffusive dynamics of oxygen and silicon is similar 
(Do I D si ~ 1.3). Upon cooling the slowdown in silicon 
dynamics is significantly stronger as compared to that in 
oxygen dynamics. This results in a faster decrease of D$i 
at low temperatures, e.g at 2800K one has Do/Dsi ~ 2.3 
|{|. This has been regarded as an evidence for a decou- 
pling dynamics. Interestingly, the temperature depen- 
dence of the average time-scale for the lifetime of Si-0 
bonds follows D (T) rather than D Si (T) @. Further- 
more, the self part of the van Hove correlation function 
for oxygen displays a clear additional peak, which corre- 
sponds to the hopping dynamics. In contrast, no clear 
hopping peak is present for the case of silicon |g ■ More- 
over, the recent study on the facilitated dynamics in liq- 
uid silica 0, IH shows that the cooperative motion for 
oxygen dynamics can be described by a string-like mo- 
tion 0] whereas the usual string-like motion is absent for 



the case of silicon atoms. 

The potential energy landscape (PEL) has been shown 
to be an indispensable tool to elucidate the nature of the 
dynamics and thermodynamics of liquids 0, 0, 0, 0, 
D_J EH EE 0, EH • In this framework the dynamics of the 
system is represented by the dynamics of a point in the 
high dimensional configurational space E3- The ther- 
modynamics of the system can be expressed in terms of 
the statistical properties of the inherent structures (IS), 
i.e. the local minima of the PEL El HI Natu - 
rally, the essential features observed in the small sub- 
system are trivially averaged out on a macroscopic level 
[2TI l22l| . Therefore, to understand the decoupling dy- 
namics in terms of the PEL properties it is important to 
analyze sufficiently small systems without relevant finite- 
size effects. 

Recent studies on the thermodynamics of this model 
have shown that the distribution of IS energies G(e) 
possesses a low-energy cutoff e c Below e c the 

configurational energy is strongly depleted as compared 
to the Gaussian distribution, which describes G(e) for 
e > e c . An important step to characterize the dynamics 
has been the introduction of metabasins (MB) [lH |23| . 
This corresponds to a coarse-graining of the configura- 
tion space. Groups of minima between which the system 
performs forward and backward jumps are regarded as 
single states, i.e MB. The total dynamics can then be 
formally described as a sequence of jumps between MB, 
each characterized by a waiting time r and an energy e 
which is chosen as the minimum IS energy of the cor- 
responding group of IS. Actually, the distribution G(e) 
of MB energies is virtually identical to that of IS such 
that the thermodynamics does not change when switch- 
ing from the IS to the MB description. It turns out for 
silica that on this level of description the oxygen dynam- 
ics of the system can be (to a good approximation) char- 
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acterized as a random- walk in configuration space with a 
temperature-independent average hopping distance [23| . 
This provides a strong link between long-range trans- 
port and microscopic dynamics. As an immediate con- 
sequence the temperature dependence of Do{T) is de- 
termined by the temperature dependence of the average 
waiting time. Furthermore, it has been shown that the 
diffusion and relaxation in liquid silica has been related 
to the defect structure j^, 0, Hf| since close to defects 
the mobility is largely enhanced. 

The goal of this paper is twofold. First, we want to 
show why the decoupling of oxygen and silicon dynam- 
ics as well as the proportionality of the oxygen diffusion 
with the inverse bond lifetime is a natural consequence 
of the nature of the dynamics of BKS-silica. It will turn 
out that it is the presence of rotational displacements 
which are the origin of the decoupling. Second, we will 
relate the occurrence of these specific displacements to 
the properties of the PEL. This will be another example 
where the reference to properties of the PEL gives new 
insight into the complex dynamics of supercooled liquids. 

The outline is as follows. In Sect. II we present some 
computational details. In Sect. Ill we show a detailed 
analysis of the nature of the MB transitions. In this way 
the relevance of rotational processes is elucidated. Based 
on these results we can unambiguously explain the de- 
coupling of the oxygen and silicon diffusion in terms of 
these processes, as shown in Sect. IV. In Sect.V the rota- 
tional processes are related to properties of the PEL. We 
close with a discussion and a summary in Sect. VI. In par- 
ticular, we show that the general physical picture of the 
relaxation processes in silica, as derived in our previous 
work |23j, is fully compatible with the results, obtained 
for the relation of the oxygen and silicon dynamics in this 
work. 



BKS-silica (in terms of the diffusive activation energy 
and the degree of non-exponentiality in the incoherent 
scattering function) in the accessible temperature range 
and the properties of tunneling systems [22l I23I l24l l2q . 
As in our previous work we use a system size of N = 99. 
To characterize the distance between two configurations 
we use the displacement of the oxygen atoms if not men- 
tioned otherwise. 

For each equilibrium run, the system is quenched via 
the standard conjugated-gradient method. The times 
for the quench have been optimized via the bisection- 
ing method to enable an efficient search for the precise 
location of the metabasin transitions |T^ | . In this way we 
obtain a time series of IS energies, which contains all rel- 
evant IS transitions which are necessary to characterize 
all MB transitions. As mentioned above MB are charac- 
terized by a waiting time r and an energy e. The energy 
of a MB is defined as the lowest IS energy in a MB (see 
[27j for a precise definition of MB). 

For later purposes we introduce two different energy 
distributions, relevant at temperature T(= l/(feg/3)). 



p(e,T) oc G(e) exp(— j3e) 



(1) 



is the Boltzmann weighted distribution of MB energies 
at temperature T j2^|. It expresses the probability that 
at a randomly given time the system is in a MB of energy 
e. Furthermore, we define 



if{e,T) 



p(e,T)(r(T)) 
(r(e,T)) 



(2) 



where (r(e, T)) is the average waiting time of a MB with 
energy e. tp(e, T) denotes the probability that after a MB 
transition the new MB has the energy e [T^ . 



II. TECHNICAL DETAILS 

We have performed molecular dynamics simulations 
of liquid silica at constant volume and energy (micro- 
canonical ensemble). The density of liquid silica is 
p = 2.30 gcm~ 3 , which is similar to that in previous 
simulations 0, ^| . The standard BKS potential @ is 
used to describe the interaction energy for liquid sil- 
ica. The Coulombic interaction is calculated with the 
Ewald-summation technique. The short range interac- 
tion is truncated and shifted to zero at the cutoff radius 
r c . To avoid an energy drift r c is set to 8.5 A, which is 
slightly larger than the one used previously |5j. In this 
work we analyze the temperature range from 2800K to 
4000K, using well-equilibrated samples. All results have 
been averaged over three independent runs in order to 
gain better statistics. 

It has been shown that for BKS-silica a system size 
of approx. 100 particles with periodic boundary condi- 
tions is sufficient to avoid relevant finite-size effects for 
the configurational entropy, the relaxation dynamics of 



III. DYNAMIC PROPERTIES OF OXYGEN 
AND SILICON ATOMS 

A. Metabasin analysis 

In the previous studies on the binary Lennard- Jones 
system (BMLJ) it has been shown that the mean square 
displacement R 2 (umb) after umb transitions between 
MB is independent of temperature except for extremely 
high temperatures for which minima of the PEL are 
not relevant due to the dominance of anharmonic ef- 
fects 01 • As a consequence the diffusion constant D(T) 
is related to the average MB waiting time (r(T)} via a 
temperature-independent jump length d, i.e. 

D(T) = d 2 /6N(r(T)). (3) 

Here we repeat this analysis for the oxygen and sil- 
icon dynamics in silica by calculating Rq(umb) and 
R-Si( n MB), respectively; see Fig^ One directly observes 
that R^j^mb) is temperature- independent whereas 
there is a significant and systematic temperature shift 
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Figure 1: The average mean square displacement of oxy- 
gen and silicon after n MB transitions, i.e. 7?o(nA/s) and 
Ft? Si {nMB)- The value of umb indicates the number of MB 
transitions. For the weak but significant temperature depen- 
dence of Rsi(riMB)/nMB in the constant regime for large 
Umb one obtains 1.3 A 2 , 1.2 A 2 , 0.9 A 2 and 0.8 A 2 for the four 
temperatures (starting from high temperatures) , respectively. 



of R 2 Si (riMB)- Actually, the temperature independence 
of Rq(umb) has been already used in our previous work 
[23| to justify the use of Eq|3] also for silica. Of course, 
the decoupling of oxygen and silicon directly follows from 
the different behavior of Ro(umb) and R 2 sii n MB)- 

On a qualitative level these observations agree with 
the previous analysis of the decoupling of oxygen and 
silicon 5j. By analyzing the time decay tb of the Si- 
O bond residual time, Horbach and Kob observed that 
tb • Do does not depend on temperature whereas tb ■ Dgi 
shows a significant temperature dependence at low tem- 
peratures, tb should be proportional to the waiting time 
of MB tmb because typical MB transitions will involve 
bond-breaking events |23. Thus the previous observa- 
tions Dq 1 oc tb but -Dgj 1 9^ tb are indeed qualitatively 
equivalent to the results in FigQ 

But how to rationalize the temperature dependence 
°f -^si( n Ms)? The most simple explanation would be 
that typical silicon displacements during MB transitions 
are smaller at lower temperatures. In this case the tem- 
perature dependence would reflect a change of the jump 
mechanism with temperature. In literature, the com- 
plicated relaxation of silicon atoms has been related to 
the possible presence of deep potential traps for silicon 
0. In this scenario, however, one would expect that 
there is also a systematic temperature dependence for, 
e.g. R 2 si{nMB — 1), which corresponds to the typical 
displacement after one MB transition. Interestingly, in 
our simulations we have not found such a temperature 
dependence; see Fig^ 



B. Rotational and translational periods 

Here we show that a different scenario is responsible 
for the observations, shown in FigQ For each time series 
of IS (MB) at different temperatures, we have observed 
some rotational processes. More precisely this means 
that during the time evolution the system is found in 
different configurations which are identical except for a 
permutation of oxygen atoms. The most simple realiza- 
tion of rotational processes are rotations of a Si04 tetra- 
hedral unit along the three-fold (C3) or two-fold (C2) axis 
but also more complex permutations have been observed. 

To elucidate this possibility in a quantitative way we 
first define rotational periods as the maximum time in- 
tervals [t\,t2\ with the property that the configurations 
at t\ and t^ are identical except for the permutation of 
oxygen atoms. Of course, during one rotational period 
several rotations may take place. To characterize the dy- 
namics during a rotational period further one can iden- 
tify the different oxygen permutations. They are denoted 
rotational states. Naturally, each rotational period con- 
tains at least two rotational states but also more may be 
present. 

In Fig|3 a sequence of IS is shown. The different rota- 
tional states which are visited during this time period are 
indicated. Note that within the rotational period the sys- 
tem may escape from a rotational state and either come 
back to the same state or end up in another permutation 
of oxygen atoms. Typically, subsequent rotational states 
are not connected by a single saddle but rather by high- 
energy IS. We note in passing that for all simulation runs 
we have not observed any permutation of silicon atoms. 
The time intervals between the rotational periods are de- 
noted translational periods; see Fig|3 Whereas the con- 
tributions from the long-time dynamics of oxygen atoms 
come from rotational and translational periods, the long- 
time dynamics of silicon atoms is exclusively related to 
the translational periods. 

On average, many IS transitions are required for the 
system to go from one rotational state to another; see 
FigI3 During this time also the silicon atoms display 
significant motion but finally go back to the original 
position. From a single-atom perspective during ro- 
tational processes the position of some oxygen atoms 
change. Thus, rotational processes contribute to the 
long-range dynamics of oxygen as opposed to that of sil- 
icon. The simulation data show that rotational periods 
become more relevant at low temperatures (see below 
for a quantification of this statement). Thus, the frac- 
tion of MB transitions which finally are irrelevant for 
the silicon dynamics increases with decreasing tempera- 
ture. From this observation it is indeed expected that 
Rsi( n MB) is smaller at lower temperatures as observed 
in Fig^ Furthermore, a much weaker temperature de- 
pendence should occur for R^^mb = 1) because on the 
single-transition scale a silicon atom hardly realizes that 
a few MB-transitions later it will go back to the origi- 
nal configuration. Stated differently, the strong forward- 
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Figure 2: The time series of IS energies at 3500K. The trans- 
lational and rotational periods are labeled by trans and rot, 
respectively. 
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Figure 3: The average mean square displacement of silicon 
RliB,Si( n ) after n MB transitions. In contrast to Fig.l only 
MB transitions within the translational periods are taken into 
account. 



backward correlations of the silicon dynamics during the 
rotational periods should show up by a stronger subdif- 
fusive behavior in R^^Umb) m agreement with Fig^ 

The hypothesis that the slowing-down of silicon is ex- 
clusively due to the presence of rotational periods can 
now easily be tested. For this purpose we have calcu- 
lated R 2 Si (n tr ans) by ignoring the MB transitions during 
the rotational periods. Thus for this analysis we only 
count MB transitions during the translational periods, 
as indicated by the index trans. The result is shown 
in Fig|3 Within statistical uncertainties one observes a 
temperature-independent behavior. This unambiguously 
shows that the decoupling of oxygen and silicon dynamics 
is exclusively due to the increasing relevance of rotational 
processes at lower temperatures. Furthermore, the non- 
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Figure 4: Temperature dependence of the rotational and 
translational spatial factors d^ rans Q (T), d^ rans ^ Si (T) and 
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for MB transitions (see definitions in text). 



diffusive part for small n tra ns, reflecting deviations from 
a pure random- walk behavior, becomes smaller and simi- 
lar in size to that of the oxygen dynamics in Fig^ From 
the large-n tr , Q „ s behavior of R 2 Si {n trans ) cx n trans the pro- 
portionality constant can be interpreted as an effective 
jump-length d 2 rans Si . It is shown in FigQ] As expected, 
no relevant temperature dependence is observed. 

Is the oxygen dynamics different in rotational as com- 
pared to translational periods? Starting with the oxygen 
dynamics during rotational periods we first determine the 
number of MB transitions n ro t,i and the squared displace- 
ment d 2 during the i-th rotational period. Formally, d\ is 
the squared Euclidian distance between the first and the 
last configuration of a rotational period. Then one can de- 
fine the effective oxygen hopping distance d 2 ot during 
rotational periods via 



d 2 

a rot,0 



n rot 



(4) 



Its temperature dependence is also shown in Fig0| 
Again, within statistical uncertainties no temperature de- 
pendence is observed. 

Formally, the total displacement of the oxygen atoms 
during the translational periods can be obtained from 
the total mean square displacement during the whole 
simulation run minus the contributions from the rota- 
tional periods, i.e. ^\ d 2 . For the number of MB tran- 
sitions ntrans in the translational periods one simply has 
ntrans = n M b — "rot- In analogy to EqQ] one can de- 
fine d 2 rans as the effective hopping distance during the 
translational periods. Again, as seen in Fig0] also this 
observable does not display any significant temperature 
dependence. 

The ratio of d 2 trans / d 2 trans Sl is close to three which 
is larger than the statistical factor of 2. Thus, the sili- 
con dynamics is indeed somewhat more hindered than the 
oxygen dynamics. This reduction in mobility, however, is 
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Figure 5: Temperature dependence of d 2 aucc rot (see definition 
in text). 

temperature-independent for the range of temperatures, 
analyzed in this work. Apart from the temperature- 
independence one observes a surprising similarity of 
d% ot o an d ^transO- This implies that, maybe in con- 
trast to expectation, the nature of the oxygen dynamics 
in translational time intervals is similar to that in rota- 
tional time intervals. 
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Figure 6: Top: The distribution of the distances between 
successive rotational states at T — 3000 K as identified during 
our MD simulations. Note the logarithmic scale. 



A second but significantly smaller peak at approx. 13 
A 2 reflects the exchange of two adjacent oxygen atoms. 
Furthermore, a very small fraction of more complex per- 
mutations with larger displacements is observed. 



IV. QUANTIFICATION OF THE DECOUPLING 
OF OXYGEN AND SILICON DYNAMICS 



C. Transitions between rotational states 

In the final part of this section we discuss some prop- 
erties of the transitions between rotational states. These 
results will shed some light on the relation between the 
PEL-properties and the actual dynamics, as will be dis- 
cussed in Sect. VI. 

First, we analyze the typical distance in the configura- 
tion space between successive rotational states. For this 
purpose we define 

.2 = I2i d l (t -\ 

u succ,rot,0 — \°) 
't j succ,rot 

where n succ ^ rot denotes the number of transition between 
different rotational states, added over all rotational peri- 
ods. Its temperature dependence is shown in FigEl Upon 
cooling d 2 succ rot increases. This means that the typical 
distance between subsequent rotational states becomes 
larger for lower temperatures, i.e. more complex oxygen 
permutations occur. As an immediate consequence on 
average more MB transitions are necessary at lower tem- 
peratures to reach a new rotational state within the same 
rotational period. 

More detailed insight into the rotational dynamics can 
be obtained by analyzing the distribution of the mean 
square displacement between two successive rotational 
states (here: T = 3000 K). This distribution is shown in 
FigEl One can clearly see that most displacements are 
around 20 A 2 which corresponds to a simple C3 rotation. 



In this section we first analyze the ratio of the oxy- 
gen diffusion constant to the silicon diffusion constant 
Do(T)/Dsi(T). Its temperature dependence, as ob- 
tained from our molecular dynamics simulations, is 
shown in Fig[7J Included are also data from the pre- 
vious MD work by Horbach et al [f| and by Saika-Voivod 
et al [l3T | . Furthermore we checked for two temperatures 
that for our system parameters the ratio does not de- 
pend on the system size within statistical accuracy. Note 
that the mild variations of this ratio for different simula- 
tions may be due to the slightly different densities or the 
different cutoff conditions for the pair-potentials. In any 
event, on a semi-quantitative level all simulations display 
a very similar increase of Do(T) / D si(T) with decreasing 
temperature. 

Using the quantities, introduced in the previous sec- 
tion, we may write 

D (T) n trans (T) ■ d 2 trans 

d r ot,o 

D Sl (T) ! 2 ■ n trans (T) • cP tranStSi ' 

trans,0 

(T) , ^rot.O 
2d tr ans,Si ^ n trans ' d trans Si 

(6) 

The factor of 2 reflects the stoichiometry of SiC>2. Using 
the results from Figgl i.e. d 2 trans / d 2 trans Sl w 3 and 

d-lot.o ~ dt r ans,o> tne ratio Dq/D S i can be rewritten as 
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Figure 7: The decoupling between oxygen and silicon dy- 
namics. Shown is the temperature dependence of the ratio 
of the oxygen diffusion constant to the silicon diffusion con- 
stant Do I D si as well as the theoretical estimations from Eq.7 
and Eq.8. Moreover, the corresponding values from literature 
[Mll^l are included. 



Do(T) 
D Sl {T) 



n ro t{T) 



(7) 



This result quantifies the general statement, formu- 
lated above, that the decoupling is exclusively related 
to the increasing relevance of rotations at lower tem- 
peratures, as expressed by the temperature depen- 
dence of n rot (T)/n t rans(T). The estimated values of 
Do{T)/Dsi{T) according to Eq[7| are also included in 
Fig El The excellent agreement with the simulated 
Do(T)/Dsi(T) constitutes a consistency check of the 
analysis of Sect. III. 



V. ENERGETIC ASPECTS OF MB 
TRANSITIONS 

MBs are mainly characterized by energy and residence 
time. In literature it has been shown that the energy 
plays an important role for the dynamics [UI27]]. Thus, 
one may expect that energy also affects the nature of the 
MB transitions in the present context. In order to check 
the energy dependence for the rotational and transla- 
tional transitions, we determine the number of MB tran- 
sitions in the rotational periods n rot (e, T) and the trans- 
lational periods n tra ns(e, T) in dependence on MB energy 
e. From Eq[7\D Sl (T)/D (T) can be written as 



DsijT) 
Do(T) 



in. 



3 Tltrar. 



S (T) + n rot (T) 
^ J de <p(e, T) ■ P tr ans(e, T) 



(8) 



ip(e, T) has been defined in Eq[3 and denotes the 
probability that a MB has the energy e. Furthermore, 
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Figure 8: Energy (a) and temperature (b) dependence of the 
probability of translational relaxations Ptrans{e.,T). 



Ptrans{e,T) is the probability that a MB of energy e is 
visited during a translational period. A small value of 
Ptrans(e,T) thus implies that after leaving a MB of en- 
ergy e the system very likely will later on end up in a ro- 
tational state. The energy and temperature dependence 
of Ptrans{e, T) is shown in Figs|3a) and (b), respectively. 
Most importantly, Pt rans (e,T) depends strongly on en- 
ergy whereas it does not show any significant tempera- 
ture dependence. In Fig|SJa) one can separate the en- 
ergy dependence of Pt rarls (e, T) into two regimes. In the 
high-energy regime one has Ptrans(e, T) m 1. This means 
that no significant rotational processes happen for high- 
energy MB. When the MB energy decreases, rotational 
processes become relevant. This results in a decrease of 
Ptrans(e,T). Interestingly, more than 90% of the tran- 
sitions involving MB with energies close to the cutoff 
energy e c contribute to the rotational dynamics. From 
the present results we have to conclude that energy not 
only determines relaxation rates as shown in the previ- 
ous studies H3,|23 but also determines the relevance of 
rotational processes. 

As a consistent check, we calculate Do{T)/ Dsi{T) via 
Eq. [S] The results are also included in Fig[7| Nat- 
urally, a good agreement with the actually observed 
D {T)/D Sl {T) is found. 
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Figure 9: The distribution tp(e, T) of MB for different temper- 
atures, fn the inset the temperature dependence of ifiom = 
/g° +le de<p(e,T) is shown. It reflects the number of low- 
energy MB close to the cutoff of the PEL. 

VI. DISCUSSION AND SUMMARY 

In previous work on silica we were able to express the 
oxygen diffusion constant via 

D (T)k J dep(e,T)/(r(e,T)}. (9) 

where the average waiting time (r(e,T)) can be written 
as 

(T(e,T)) = T (e)exp(/3V(e)) (10) 

with an energy-dependent attempt frequency l/ro(e) and 
an effective barrier height V(e). 

This set of equations implies that from the three 
temperature-independent quantities G(e),To(e) and V(e) 
the thermodynamics as well as the temperature depen- 
dence of the oxygen dynamics (in terms of the diffusion 
constant) can be predicted. In particular V(e) can be ra- 
tionalized in terms of the local topology of the PEL • 
This implies a basically complete understanding of the 
oxygen diffusion in terms of properties of the PEL. In this 
work we have identified a further (nearly) temperature- 
independent function -Pt rtms (e) which allows one to pre- 
dict the enhanced slowing-down of silicon diffusion as 
compared to oxygen diffusion. 

In the following we will argue that the results, obtained 
so far, are compatible with the general physical scenario 
of relaxation processes in silica, as formulated in [23j |. 
As a main result it has been shown that there exists a 
crossover-energy e cross = -1905 eV which separates two 
dynamical regimes. For e > e cross the escape dynamics is 
non-activated whereas for significantly lower e the system 
first has to climb up in the PEL to IS with energies close 
to e cross until it can effectively escape from this region 
of configuration space. On a qualitative level this obser- 
vation can be related to the strong energy dependence of 
the density of MB G(e) or <p(e,T), respectively. In the 
considered range of temperatures much less than 1% of all 



configurations possess an energy less than e c + 1 eV; see 
Fig|5] In contrast, nearly every second MB has an energy 
less than e cross . As a consequence for MB with e > e cross 
the PEL supplies many adjacent configurations of similar 
energy with connecting saddles of at most 1 eV which the 
system may easily reach, giving rise to non-activated dy- 
namics whereas in the low-energy limit this is no longer 
possible due to the lack of such configurations. 

For rationalizing the features, observed in this work, 
we discuss in qualitative terms the distribution of MB 
in configuration space relative to a given MB (denoted 
central MB) with energy eo, visited at a given tempera- 
ture T. We introduce the function ip a u(e, d 2 |e ) as the 
probability that a randomly selected MB is a distance d 2 
away from the central MB and has an energy e. This 
probability distribution has two contributions 

(p a ii{e,d 2 \e ) = p ro t(e, d 2 \e ) + ip rest (e, d 2 \e ). (11) 

The first contribution takes into account configurations 
which result from oxygen permutations of the central 
MB. This distribution can be written as ip ro t(e,d 2 \eo) oc 
G rot {d 2 )S(e — eo) where G rot (d 2 ) is the distribution of 
distances, when taking into account all possible 66! per- 
mutations. We have calculated G ro t(d 2 ) for one represen- 
tative low-energy configuration and have exchanged the 
indices of the oxygen atoms via Monte-Carlo techniques. 
To determine the distribution over many orders of magni- 
tude we have used appropriately adjusted flat-histogram 
techniques 29] . The resulting distribution G ro t(d 2 ) is 
shown in FigEH One can clearly see a few peaks for 
smaller distances and a transition to a continuous be- 
havior for larger distances. Naturally, the <i 2 -values of 
the peaks agree with those observed during the actual 
transitions between rotational states. For d 2 > 100 A 2 
this distribution no longer shows any wiggles. This re- 
flects the fact that a transition with a large value of d 2 
can be achieved by a multitude of different oxygen per- 
mutations. In contrast, for small d 2 discrete distances 
can be identified which belong to specific types of reori- 
entations (like C3 rotations). Comparison with the actu- 
ally observed transitions in Fig|S|shows that in particular 
more complex permutations, related to larger values of 
d 2 , are strongly suppressed. Interestingly, also exchange 
processes of adjacent oxygen atoms, corresponding to the 
first peak, are significantly suppressed. 

The second contribution ip rest (e,d 2 \eo) is related to 
the remaining MB which have a different structure as 
the central MB. For sufficiently large d 2 one may natu- 
rally write iprest(e,d 2 \e ) tx G rest (d 2 )ip(e,T) with (p(e, T) 
defined in Eq|21 Qualitatively, this means that for re- 
gions sufficiently far away from the central MB the dis- 
tribution of MB is not influenced by the presence of a 
MB with e — eo at the origin. One expects that for 
large d 2 the distribution G rest (d 2 ) increases much faster 
than G ro t(d 2 ), i.e. with a larger exponent of d 2 , be- 
cause G res t(d 2 ) contains the exponentially larger num- 
ber of different IS as well their permutations (except for 
those taken into account by G ro t(d 2 )). Thus, for large d 2 
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Figure 10: The distribution G ro t(d 2 ) of all distances between 
configurations with oxygen permutations, as obtained from 
Monte-Carlo simulations. 



one expects G rot (d 2 ) <C G rest (d 2 ). This dramatic differ- 
ence may disappear for small d 2 because the ratio of both 
functions becomes smaller for smaller d 2 . Furthermore, 
in contrast to the oxygen permutations the possible tran- 
sitions to adjacent MB are more continuously distributed 
with respect to their distances d 2 . As a consequence the 
continuous regime of G res t(d 2 ) will probably extend to 
much smaller values of d 2 than observed for G ro t(d 2 ), i.e. 
approx. 100 A 2 . Thus, one may speculate that for small 
d 2 (e.g. d 2 « 20 A 2 ) G rot {d 2 ) is of the same order as 
G res t(d 2 ). Of course, for a more quantitative character- 
ization one directly has to determine G res t(d 2 ), which, 
however, is a formidable task. 

We note in passing that the factorization of 
ip res t(e,d 2 \eo) tx G res t(d 2 )(p(e,T) will break down for 
d 2 — > 0, as known from the study of the tunneling sys- 
tems [30| . There it turned out that for very nearby IS 
(d 2 < 1 A 2 ) the energy difference does not depend on eo 
which implies that ip rest (e, d 2 |eo) is a function of e — eo. 

With the general concepts, introduced above, the fol- 
lowing observations have to be rationalized: (i) The pres- 
ence of rotations with small values of d 2 ; see Fig. 
but the strong suppression of rotations with larger d 2 . 
(ii) The additional suppression of particle exchange with 
d 2 w 13 A 2 , (in) The relevance of rotations only for 
e < e cross and the increase of rotations for decreasing 
energy, (iv) The increase of d 2 ucc rot with decreas- 
ing temperature, (v) The temperature-independence of 
Ptrans(e)- (vi) The similarity of d 2 0t and d 2 trans . 

(i) For large d 2 the function G res t{d 2 ) will be orders of 
magnitude larger than G ro t(d 2 ). Thus for purely statisti- 
cal reasons hardly any permutations with large d 2 are ob- 
served; see FigEl In this context it is surprising that nev- 
ertheless a very small number of permutations with d 2 up 
to 80 A 2 are observed which from this purely statistical 
consideration should be extremely unlikely. Maybe there 
are specific pathways on the PEL which enhance the 
probability of specific permutations, (ii) For oxygen ex- 



change processes massive reorientations are necessary in 
contrast to, e.g., effective C3 rotations. To perform an ex- 
change process the system would have to reach regions of 
the configuration space which are far away from the cen- 
tral MB. As a consequence, for purely statistical reasons 
the probability to complete the exchange process is very 
small, (iii) As discussed above for small d 2 the function 
G r ot(d 2 ) may be of the same order as G res t(d 2 ). Further- 
more, Lp res t{e, d 2 \eo) (x tp(e,T) is very small for e w e c , 
see FigEl implying (p a ii{e , d 2 \e ) « G ro t(d 2 )¥w(e , d 2 ) 
for eo ~ e c . Thus, for statistical reasons the relevance 
of rotational states is significantly enhanced for very low 
energies. In contrast, for eo > e cross the system pos- 
sesses a large number of adjacent IS with similar energy 
such that again for statistical reasons the probability to 
reach one of the rotational states is very small, (iv) As 
a further consequence the range of distances for which 
the probability to find a rotational state is non-negligible 
increases with decreasing energy. Thus, in agreement 
with observation one would expect an increase of the av- 
erage distance between successive rotational states with 
decreasing energy and thus with decreasing temperature, 
(v) As discussed above for the transition from one low- 
energy configuration (eo ~ e c ) to a new low-energy con- 
figuration the system first has to climb up the PEL to 
energies close to e cross . This trajectory basically forms 
an energy- upwards path without major intermediate bar- 
riers |23| . In analogy, the subsequent path to a new low- 
energy configuration is not expected to possess significant 
barriers, in similarity to the funnel-like picture of protein 
folding. Stated differently, the search for the subsequent 
low-energy configuration is basically entropy-driven so 
that Ptrans, indeed, should not significantly depend on 
temperature, (vi) Whether or not the system ends up 
in a rotational state is mainly a question of the statis- 
tical availability of these states. As a consequence one 
would expect that the nature of the oxygen dynamics in 
translational and rotational periods should not differ too 
much. This is consistent with the observed similarity of 

W 2 and rl 2 

u rot,0 AnLL u trans,0- 

According to Eq[S] the fraction of the rotational and 
translational dynamics is related to the thermodynamic 
weighting of the MB population, ip(e,T). Therefore, the 
temperature dependence of n tra ns (T) / Um b (T) presented 
above is just a consequence of the temperature depen- 
dence of <p(e, T). It is interesting to estimate the value of 
Do/Dsi at the glass transition temperature T g (1450K). 
From the density of MBs G(e) and waiting time (r(e, T)) 
one can predict ip(e,T) at T g . By using Eq0 the esti- 
mated value for Do/Dsi at T g is of the order of 10 1 , 
which is one order of magnitude smaller than the exper- 
imental value of « 10 2 2,0|. This discrepancy may be 
resolved in different ways: (i) At lower temperatures fur- 
ther decoupling mechanisms set in. (ii) The temperature- 
independence of Ptrans no longer holds for temperatures 
closer to the calorimetric glass transition, (iii) The limits 
of the BKS-potential to predict the details of the oxygen 
and silicon dynamics in silica are reached. 
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The present PEL approach to characterize the dynam- 
ics of oxygen and silicon allows one to obtain a relatively 
simple physical picture about the oxygen dynamics as 
well as the silicon dynamics. The essential key ingredi- 
ent to understand the onset of the decoupling of oxygen 
and silicon diffusion is the presence of rotational dynam- 
ics which only contributes to the oxygen dynamics and 
which can be rationalized in simple statistical terms in 



the framework of the PEL approach. Excluding the ro- 
tational periods, the difference in activation energies of 
silicon and oxygen disappears in the considered range of 
temperatures. 

We would like to acknowledge discussions with M. Vo- 
gel on this subject and the NRW Graduate School of 
Chemistry and the DFG (SFB 458) for the support of 
this work. 
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